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Abstract 



One of the most important topics in geophysics is to study convection in a sea. Based on the 
algebraic characteristics of the equations of dynamic convection in a sea, we introduce various schemes 
with multiple parameter functions to solve these equations and obtain families of new explicit exact 
solutions with multiple parameter functions. Moreover, symmetry transformations are used to simplify 
our arguments. 

1 Introduction 

Both the atmospheric and oceanic flows are influenced by the rotation of the earth. In fact, the 
fast rotation and small aspect ratio are two main characteristics of the large scale atmospheric 
and oceanic flows. The small aspect ratio characteristic leads to the primitive equations, and 
the fast rotation leads to the quasi-geostropic equations (cf. [2], [7], [8], [10]). A main objective 
in climate dynamics and in geophysical fluid dynamics is to understand and predict the periodic, 
quasi-periodic, aperiodic, and fully turbulent characteristics of the large scale atmospheric and 
oceanic flows (e.g., cf. [4], [6]). 

The general model of atmospheric and oceanic flows is very complicated. Various simplified 
models had been established and studied. For instance, Boussinesq equations are simpler models 
in atmospheric sciences (e.g., cf. [9]). Chae [1] proved the global regularity, and Hou and Li [3] 
obtained the well-posedness of the two-dimensional equations. Hsia, Ma and Wang [4] studied 
the bifurcation and periodic solutions of the three-dimensional equations. 

The following equations in geophysics 



U:j: +Vy + Wz = 0, 



P 



Pz 




Pt + UPX + VPy + WPz 
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Ut + UUx + VUy + WUz + V 



2000 Mathematical Subject Classification. Primary 35Q35, 35C05; Secondary 35L60. 
^Research supported by China NSF 10871193 



Vt + UVx + VVy + WVz — U = Py, (1-4) 

are used to describe the dynamic convection in a sea, where u, v and w are components of 
velocity vector of relative motion of fluid in Cartesian coordinates {x,y,z), p = p{x,y,z,t) is 
the density of fluid and p is the pressure (e.g., cf. Page 203 in [5]). Ovsiannikov determined the 
Lie point symmetries of the above equations and found two very special solutions (cf. [5]). 

In [11], we used the stable range of nonlinear term to solve the equation of nonstationary 
transonic gas flow. Moreover, we [12] solved the three-dimensional Navior-Stokes equations 
by asymmetric techniques and moving frames. Based on the algebraic characteristics of the 
equations (1.1)-(1.4) of dynamic convection in a sea, we introduce various schemes with multiple 
parameter functions to solve these equations and obtain families of new explicit exact solutions 
with multiple parameter functions. Moreover, symmetry transformations are used to simplify 
our arguments. By specifying these parameter functions, one can obtain the solutions of certain 
practical models. 

For convenience, we always assume that all the involved partial derivatives of related func- 
tions always exist and we can change orders of taking partial derivatives. We also use prime 
' to denote the derivative of any one-variable function. We will use the following symmetry 
transformations T1-T4 due to Ovsiannikov ( cf. Page 204 of [5]) of the equation (1.1)-(1.4) to 



simplify our solutions: 

Ti{u) = u{t, X + a,y, z + ol' x — a'y) — a', Ti{v) = v{t, x + a,y, z + a''x — a'y), (1.5) 

Ti{w) = w{t, X + a,y, z + ol' x — a'y) — ol' u + olv — a'' x + a''y, (1-6) 

T2{u) = u{t, x,y + a, z + a'x + a''y), ?2(f ) = v{t, x,y + a, z + a'x + a''y) — a', (1.7) 

T2{w) = w{t, x,y + a, z + a'x + a''y) — a'u — a''v — a''x — a'''y, (1-8) 

Ti{p)=p{t,x + a,y,z + a"x-a'y), T2{p) = p{t,x,y + a, z + a'x + a"y), (1.9) 

T3{u) = u{t,x,y,z + a), T^^v) = v{t, x,y, z + a), (1-10) 

T3{w)=w{t,x,y,z + a)-a', T^(jp) = p{t,x,y,z + a), (1.11) 

Ti{p)=p + a, Ti{F)=F ioT F = u,v,w, (1.12) 



where a is an arbitrary function of t. The above transformations transform one solution of the 
equations (1.1)-(1.4) into another solution. Applying the above transformations to any solution 
found in this paper will yield another solution with four extra parameter functions. 

In Section 2, we use a new variable of moving line to solve the equations (1.1)-(1.4). An 
approach of using the product of cylindrical invariant function with z is introduced in Section 
3. In Section 4, we reduce the three-dimensional (spaeial) equations (1.1)-(1.4) into a two- 
dimensional problem and then solve it with three different ansatzes. 
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2 Moving-Line Approach 

Let a and (3 be given functions of t. Denote 

w = a'x + P'y + z. (2.1) 
Suppose that /, g, h are functions in t, x, y, z that are hnear in x, y, z such that 

fx + 9y + h, = 0. (2.2) 

We assume 

u = (t){t,w) + f, v = il;{t,w)+g, (2.3) 

w = h- a'<l){t,w) - I3'i}{t,w), p = C{t,Tu), (2.4) 

where (p,ip,C are two-variable functions. Note that the first equation in (1.1) naturally holds 
and p = Pz = Cvj the second equation in (1.1). Moreover, (1.2)-(1.4) become 

C^t + C^Ua"x + P"y + a'f + P'g + h) = 0, (2.5) 

ft + 9 + ffx + gfy + hfz + a' + ^t + {fx - afz)(l> + {fy - P' + 1)V' 
+(j)Ua"x + p"y + a'f + 0g + /i) = 0, (2.6) 

gt - f + fgx + ggy + %^ + /3' + + {gx - a'g^ - 1)0 + {gy - P'gz)^ 

+i^^{a"x + 0'y + a'f + 0g + /i) = 0. (2.7) 

In order to solve the above system of partial differential equations, we assume 

a"x + 15" y + a'f + p'g + h = -j'w = --/'{a'x + p'y + z) (2.8) 

for some function 7 of t, and 

ft + 9 + ffx + 9fy + hfz + a' = 0, (2.9) 

9t-f + fgx + ggy + hgz + P' = o. (2.10) 

Then (2.5)-(2.7) become 

C^t-7'^C^^ = 0, (2.11) 

<^t + {fx - a'fz)<t> + {fy - 0fz + l)^ - ivJ<l>^ = 0, (2.12) 

V-t + {gx - a'gz - 1)<^ + {gy - P'gz)^ - I'rui/j^ = 0. (2.13) 
According to (2.8), 

h = -a"x-p"y-a'f -f3'g--f'w. (2.14) 
Substituting the above equation into (2.9) and (2.10), we have: 

ft + f{fx - a'f,) + g{fy - (3'fz + 1) - fz{a"x + p"y + -f'w) + a' = 0, (2.15) 



gt + f{gx - a'g, - 1) + g{gy - f3'g,) - g,{a"x + f3"y + ^w) + (3' = 0. (2.16) 
Our linearity assumption implies that 



A 



fx - a'f, fy - p'f, + 1 
9x - a'gz - 1 9y- P'gz 



is a matrix function of t. In order to solve the system (2.12) and (2.13), and the system 
and (2.16), we need the commutativity of A with dA/dt. For simplicity, we assume 

fy - 13' f, + 1 = 5a. - a'g, -1 = 0. 

So 

fy = P'fz - 1, 9x= a'gz + 1- 
Moreover, (2.15) and (2.16) become 

ft + fifx - oi'f,) - f,{a"x + 0'y + iw) + a' = 0, 

gt + g{gy - P'gz) - gz{a"x + p"y + j'w) + p' = o. 

Write 

/ = aix + {P'a2 — l)y + a2Z + a^, 

g = {a' (52 + l)x + (5xy + p2Z + Pa 

by our linearity assumption and (2.19), where ai and (3j are functions of t. 

Now (2.20) is equivalent to the following system of ordinary differential equations: 

a'l + ai(ai — ol 0.2) — a2{ot'' + j'ot') = 0, 

{(3'a2y + {P'a2 - l)(ai - a'aa) - a2{P" + j'P') = 0, 
(^2 + C(2{ai — aa2 — 7') = 0, 
+ Q!3(q!i — a' 0:2) + a' = 0. 
Observe that (2.25) - /?' x (2.26) becomes 

—ai + a'a2 = 0. 

So (2.26) becomes 

a2 — 7'a2 = a2 = hie* , hi G M. 

According to (2.28), 

With the data (2.29) and (2.30), (2.24) naturally holds. By (2.27), we take 

0:3 = —a. 
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Note that (2.21) is equivalent to the following system of ordinary differential equations: 

a'P'^ + {a' (32 + l)(/?i - 13' (32) - a'(32j' = 0, (2.32) 

(3[ + /3i(/3i - /3'/32) - P2{P" + (3'i) = 0, (2.33) 

(3'2 + (32{(3i-(3'(32-i) = Q, (2.34) 

(3'^ + (3m-(3'(32)-(3' = Q. (2.35) 

Similarly, we have: 

(3i = h2P'e\ (32 = h2e\ (3^ = (3 (2.36) 

with 62 e M. Moreover, (2.2) gives 7' = by (2.14), (2.28) and (2.36). We take 7 = 0. 
Therefore, (j) = '^{w) and ^ip = i{vj) by (2.12) and (2.13) for some one- variable functions S and 
i. Furthermore, we take = a{w) by (2.11) for another one-variable function a. In summary, 
we have: 

Theorem 2.1. Let a, (3 be functions of t and let 61,62 G Suppose that Sj, l and a 
are arbitrary one-variable functions. The following is a solution of the equations (1.1)-(1.4) of 
dynamic convection in a sea: 

u = bia'x + {bi(3' - l)y + biz-a + ^{a'x + (3'y + z), (2.37) 

V = (62a' + l)x + 62/3'y + b2Z + P + i{a'x + (3'y + z), (2.38) 

w = -{a" + bi{a'f + (62a' + l)(3')x - {(3" + a'{bi(3' - 1) + b2{(3'f)y - {bia' + b2(3')z 

+aa' - (3P' - a'^{a'x + /3'y + z) - (3'L{a'x + (3'y + z) , (2.39) 

p = a{a'x + (3'y + z), p = a'(a'x + (3'y + z). (2.40) 

We remark that we have tried some other forms of the matrix A in (2.17) such that A and 
dA/dt commute, but we have failed to get new solutions. 

3 Approach of Cylindrical Product 

Let (7 be a fixed one-variable function and set 

w = a{x'^ + y'^)z. (3.1) 
Suppose that / and g are functions in t, x, z that are linear homogeneous in y and 

h = ^-z{f, + gy), (3.2) 

where 7 is a function of t. Assume 

u = f + yi(j(t,w), V = g — xi(^{t,w), w = h, p = (j){t,w) (3-3) 



5 



where and (f) are two-variable functions. Note 

ut = ft + yipt, Ux = fx + '^xyza'tp^, 
Uy = fy + ip + 2y'^za'tp^, Uz = fz + yaip^, 
vt = gt- xipt, Vx = gx-ip- 2x^za''^^, 
Vy = Qy- 2xyza'ij;^, = gz - xai/j^. 

Hence (1.3) becomes 

V4 + uux + vuy + wuz + v = ft + y'ipt + if + y^){fx + 2xyza'il)^) 

+{g - xi'){fy + 1 + -0 + 2y^za'il)^) + yahip^ 
= ft + ffx + 9{l + fy) + x{gx-fy-l)i^-xi^^ 

2xza' 

+y[i^t + {fx + gy)ip + {2{xf + yg)a'z + ha)^;^ ' 



^roj — 



and (1.4) gives 

vt + uvx + vvy + wvz -u = gt- xtpt + (/ + yip){gx - - 2x'^za' 

+{g - xil)){gy - 2xyza''il}^) - xahip^ 

= gt + f{gx - 1) + 99y - y(i + fy- gx)'4' - yV"^ 

2yza' 

-x[i^t + Ux + gy)i^ + {2{xf + yg)a'z + ha)il)^] = —. 

In order to solve the above system of differential equations, we assume 

f = a'x-^, g = ^ + a'y, a(x^ + y^) = o ! 2 

for some function a of t. According to (3.2), 

h = -- 2a' z. 
a 

Now (3.8) becomes 

(a" + {a'f + - tp'^)x + y[i)t + 2aV + (7 - ^o'w)'4)^] = 2xw 

and (3.9) yields 

{a" + {a' f + - V'^)y - x[iJt + 2aV + (7 - 4a'tu)V'^] = 2yw. 
The above system is equivalent to 

a" + {a'f + 4-1 - = 2ro, 
i^t + 2a i) + (7 - 4a'ro)V'^ = 0. 



By (3.14), we take 

iIj = ^Ja" + {a'Y + 4-1 - 2ro, (3.16) 
due to the skew-symmetry of {u,x) and {v,y). Substituting (3.16) into (3.15), we get 

a'" + 2a'a" + Aa'{a" + {a'f + 4"^ - 2w) - 2(7 - Aa'w) = 0, (3.17) 

equivalently, 

j = 2{a'f + 3a'a" + ^-^^. (3.18) 
According to the second equation in (1.1), we have p = crcj).^. Note 

Pt = cr4'wt, Px = 2xcr'(^^ + OJ^^^), (3.19) 

Py = '^y'^'i'P^ + '^(l>vjw), Pz = o-'^(t^ww (3.20) 

So (1.2) becomes 

- 2aVCT + (7 - 4a'ro)0^j,^jj = 0. (3-21) 
Modulo T4 in (1.12), the above equation is equivalent to: 

(/)t + 2a'(/)+(7-4a'tu)(/)^ = 0. (3.22) 

Set 

^ = e^^V, ^ = e2°(/). (3.23) 
Then (3.15) and (3.22) are equivalent to the equations: 

^t + {l- ^a'^)ip^ = 0, 4>t + {l- ^ot'w)^^ = 0, (3.24) 
respectively. So we have the solution 

^ = ^(f) = e-2°9 {e^'^^a" + {a'Y + 4-i - 2w^ 
for some one- variable function 55. Thus we have: 



(3.25) 



Theorem 3.1. Let a he any function oft and let ^ he arbitrary one-variable function. The 
following is a solution of the equations (1.1)-(1.4) of dynamic convection in a sea: 



V I 2z 

u = a'x-'- + yJa" + {a'Y + --^—^, (3.26) 



v = a'y + l- x^a" + {a'f + \ - (3.27) 
w={ 2{a'f + 3a'a" + ^^-^ ) [x" + y^) - 2a' z, (3.28) 





g2a 


9' 





1 2z 



4 x"^ + y"^ 



(x2-Ky2)^Q,/ + Q,2+ 1 



2z 



(3.29) 
(3.30) 



4 x-^+y-^ 
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4 Dimensional Reduction 

Suppose that u, v, ( and 77 are functions in t, x, y. Assume 

w = C,- {ux + Vy)z, p = z + r], p=l- (4.1) 
Then the equations (1.1)- (1.4) are equivalent to the following two-dimensional problem: 

ut + uux + vuy + v = -rjx, (4.2) 
vt + uvx + vvy - u = -rjy. (4.3) 

The compatibility rixy = rjyx gives 

{Uy - Vx)t + u{Uy - Vx)x + v{Uy - Vx)y + {Ux + Vy){Uy " ^o; + l) = 0. (4.4) 

Suppose that is a function in i, x, y such that 

^xx +^yy = ^ (4.5) 
(so ■!? is a time-dependent harmonic function). We assume 

U = ^xx, V = l^xy (4.6) 

Then (4.4) naturally holds. Indeed, 

Ut + UUx +VUy + V= {l^xt + 2~\€x + €y) + ^^v) x ^ (4-7) 

Vt + uVx + vvy-u= {^xt + + €y) + ^y)y ■ (4.8) 

By (4.2) and (4.3), we take 

V=-^xt-^y-l{€x+^ly)- (4.8) 

Hence we have the following easy result: 

Proposition 4.1. Let and ^ be functions in t, x, y such that (4-5) holds. The following is 
a solution of the equations (1.1)-(1.4) of dynamic convection in a sea: 

U = T^xx, V = '&xy, W = C, (4.9) 

p=l, p = ^_^^^_^^_l(^2^ + ^2^). (4.10) 

The above approach is the well-known rotation-free approach. We are more interested in the 
approaches that the rotation may not be zero. Let / and g be functions in t, x, y that are linear 
in x,y. Denote 

w = x^ + y^. (4.11) 
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Consider 

u = f + y4>{t,'aj), V = g - x(t){t,w), (4.12) 
where ^ is a two- variable function to be determined. Then 

Ux = fx + "^xycf)^, Uy = fy + (j) + 2y^(j)^, (4.13) 

Vx= 9x-(t>- 2x^(?!)c7, uy = Qy- 2xy(f)^. (4.14) 

Thus 

Ux+Vy = fx+9y, Uy-Vx = fy-gx+2{'uj^).^. (4.15) 
For simplicity, we assume 

a'x y X a'y 

J = 9 = — -— (4.16) 

•' 2q; 2' ^ 2 2a ^ ' 

for some functions a and (3 of t. Then (4.4) becomes 

a' a' 

{zU(p)^t ^37(^^7^)^^ {w(p)^ = 0. (4.17) 

a a 



Hence 

7 + S(aro) 

9 = 

w 

for some function 7 of t and one- variable function $5. 
Now (4.12), (4.16) and (4.18) imply 



Moreover, (4.2) and (4.3) yield 



4a2 4 



w 



Thus 

J\2 _ o^,^/' 



(4.18) 



n = -^-^ + il±^^^, (4.19) 
2q! 2 w 

x_a^_ (7 + Q{aw))x 

2 2a w ' ^ ' ' 



ia'f-2aa" 1\ 7'y ,0 

^ ^ +- x + ^-a;</)2 = -%, (4.21) 



{a'f-2aa" l\ 7'x . , ^ 



"=2/ H 4.^ +^).^ + yarctan|. (4.23) 



Theorem 4.2. Lei a, 7 he any functions oft. Suppose that Q is an arbitrary one-variable 
function and C, is any function int,x,y. The following is a solution of the equations (1.1)-(1.4) 
of dynamic convection in a sea: 

y {-f + Q{{x'^ + y^)a))y 
2a 2 x'^ + y'^ ' ^ ' ^ 

X _ a'y _ (7 + a((a;2 + y2)a))x 
"^"2 2a + ' ^4-^^^ 
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ol 

w = -z + C„ p=l, (4.26) 

P = z + -j + 4 j (x^ + y') + y arctan ^ (4.27) 

with w = + . 
Next we assume 

u = e{t,x), V = (l){t,x) + i;{t,x)y, (4.28) 
where e, (f) and i/j are functions in t, x to be determined. Substituting (4.28) into (4.4), we get 

(ptx + i^txv + £{<t>xx + ipxxy) + + ^y)'^x + {^x + ip){4'x + ipxy - i) = o, (4.29) 

equivalently, 

{(l>t + e(l>^ + 4>4^ - e)^ - i/j = 0, (4.30) 

(V't + £Vx + V'')x = 0. (4.31) 



For simplicity, we take 

a function of t. 
Denote 

Then (4.30) becomes 



tp = -a', (4.32) 



= (t) + x. (4.33) 



i^t + e4 - a'^)x = 0. (4.34) 
To solve the above equation, we assume 

^ = ^-#4 (4.35) 

for some functions (3 of t, and of t and x. We have the following solution of (4.34): 

^ = e"9(i?) ^4> = e"9(7?) +x^v = e"9(??) +x-a'y (4.36) 
for another one-variable function Moreover, 

£ = ^ir ill U 37) 



Note 



(/3e-")^ /3e-"(i?^t9^(^)+^t^.Q^^(^)) - ^t^xt , 
ut + uu^ + vuy + v = ^^^ ^^^^f^, ^ ay 



vt + uvx + Wj/ — u = {{a')^ — q;'')j/ + f3 — a'x. (4.39) 
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By (4.2) and (4.3), 

Wx,-ft + ^ L±M -,{J^^-^J . (4.40) 

Theorem 4.3. Let a, [3 he functions of t and let be a one-variable function. Suppose 
that -& and ( are functions in t,x,y. The following is a solution of the equations (1.1)-(1.4) of 
dynamic convection in a sea: 

'''Swi-h -''"■^m+^-'v. (4.41) 

+ ^ Mi-— --i . (4.43) 



2 V^x5J'(^) 



Finally, we suppose that a, (3 are functions of t and /, g are functions oi t,x,y that are linear 
homogeneous in x and y. Denote = ax + Py. Assume 

u = f + P(p{t,w), V = g - a4>{t,w). (4.44) 

Then 

Uy-Vx = fy- gx + {0i'^ + 0^)4>VJ, + Vy = fx + gy (4.45) 

Now (4.4) becomes 

fyt - 9xt + (a^ + ^^)'<^^ + (a^ + /32)(<^^t + (a'x + + a/ + /Jff)'^.^^) 

+ gy){fy -gx + l + {a^ + P^)(t>^) = 0. (4.46) 

In order to solve the above equation, we assume 

gx = 'P, fy = ^-'^, (4.47) 

a'x + (3'y + af + Pg = (4.48) 
for some function (p of t. The equation (4.48) is equivalent to: 

a' + af, + ^P = 0^f, = (4.49) 
P' + f3gy + a{^-l) = 0^gy = J' + - ^) _ (4.5O) 
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Now (4.46) becomes 



+ + - ) = 0. (4.51) 



(4.52) 



Thus we have the following solution: 

= 4±4e/(«^"'(^-i)+-"'^^)'^*9'(ti7), 
where 9 is an arbitrary one-variable function. Note 

2(a')^ + ('^/3)^ + 3a'/3(/7 - a((^/?)' - aa" , 2 



+ ( — ^ — + ^^]x 

\ a (3 J 

/3((y - l)ay + (3(3" - 2{(3'f - ((y - l)af - 3a/3^(y - 1) 

By (4.2) and (4.3), 



( (3{{^ - l)ay + (3(3" - 2(/302 - {{^ - l)aY - ^a(3'{if - 1) 
2 V ^2 
x2 /2(a')^ + (¥'/3)2 + 3a'^</^ - a((^/3)' - aa" ,{^- !)("' + </^/?) 



+ ip^ + [- — - ^' 



/ \ a^ I a 

(^(/?' + a((^-l)), , a + /? 



^]., + ^ (1 - ^5^) e/(«^-(-^)--^.)-^(.). (4.55) 



Theorem 4.4. Let a, /3, 99 6e functions of t and let Q be a one-variable function. Suppose 
that ( is functions in t,x,y. The following is a solution of the equations (1.1)-(1.4) of dynamic 
convection in a sea: 

u = {<p-l)y- (^:±M^ + e/("^-^(^-^)+-^^^M*9'(ax + Py), (4.56) 

a a^ + 

«^=^^^ + ^ J^ + C, P=l, (4.58) 



, y^ f P{{^ - l)ay + PP" - 2{P'f - ((y - l)a)^ - 3a/?^(y - 1) 2\ 
P = ^ + y ^ ^^-1) J 

x2 ^ 2(a0^ + (y/3)^ + Za'P^ - a{<pPy - aa" ^ 2 V , {<P - l){a' + <pP) , 
+^^^^^±1^^] + (1 - ^5^) e/("/^-^(--)+-^/^-)'^*G(ax + Pyy (4.59) 
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